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Abstract 

QCD corrections to e^e~ WW — >• qiqiq2<l2 are computed and pre- 
sented in a form which allows to impose realistic cuts on the structure of 
the observed events. QCD radiation substantially modifies the jet-jet in- 
variant mass distributions from which the value of the W mass will be 
extracted. When the range of allowed jet-jet masses is restricted, as it has 
been proposed in order to suppress non-iyVT backgrounds, the total cross 
section is also deeply affected. If all events are forced to four jets, combining 
the two partons with smallest invariant mass, and the reconstructed masses 
are required to satisfy jM^ — My/\ < 10 GeV i=l,2 the lowest order cross 
section can be reduced by more than 40%. 
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Introduction 



One of the main goals of Lep II is the measurement of the W mass with high accuracy, 
possibly of the order of 50 MeV or less 0. Two methods have been singled out as the 
most promising . The first one is based on the rapid increase of the total cross section 
at threshold and has been recently studied in detail in ref. where it has been shown 
that the optimal collider energy is about 161 GeV. The second method relies on the 
direct reconstruction of the mass from the hadronic decay products of the W using the 
decay channels 

W^W~ qiqiq2q2 (1) 

W+W- qqiu (2) 

where i = e, /i. The tree-level relevant branching ratios are given in table I. 

The threshold method will be the first one to be exploited when Lep II is turned 
on. In the meantime the number of superconducting cavities will be increased until the 
center-of-mass energy can be pushed to about 175 GeV where the much larger cross- 
section will allow Mw to be extracted by direct reconstruction. Preliminary studies 
indicate that the direct measurement will provide a more precise determination of the 
W mass than the threshold method if all decay channels are combined, while the two 
methods are expected to have comparable accuracies if only the semileptonic decays are 
used for direct reconstruction. The possibility of a further energy upgrade to about 190 
GeV is presently under consideration. This higher energy would greatly improve the 
prospects for the discovery of new particles, in particular of the Higgs boson, without 
compromising the precision expected for the measurement of M\y. 

In the case of semileptonic decays, apart from all the usual difficulties related to 
measuring jet energies and directions, the main uncertainty in the mass measurement 
derives from the reconstruction of the unobserved neutrino. These uncertainties can 
be partially eliminated by using energy-momentum conservation and possibly the ap- 
proximate equality of the and W~ masses. In the case of hadronic decays, where 
in principle all decay products can be detected, there are additional ambiguities which 
stem from the fact that two decays occur in the same event. Even in the simplest case, 
in which only four jets are detected, there are three possible ways of pairing the jets. 
If more jets are present the combinatorics becomes rapidly quite complicated. 

It has been proposed to force all events to four jets simply iterating a given recon- 
struction algorithm until only four clusters are left. It is obvious that such a scheme 
is not without perils. A jet from, say , the decay of the can be closer to a jet 
from the decay of the W~ than to any other jet from the positively charged W. If the 
stray jet carries large energy, then the reconstructed masses will be quite far from the 
W mass which is already at present known with an error smaller than 200 MeV 0. 
One could imagine that such events could be discarded, at the same time providing 
a mean of suppressing non-W^IV background. However, due to the intrinsic width of 
the W and to all experimental uncertainties, it is impossible to impose very stringent 
cuts on the difference between the measured jet-jet masses and the true mass without 

^For a general introduction to the Standard Model predictions for VF-pair production in e+e~ 
collisions see 



substantially reducing the event rate. A precise estimate of the effect of these mass 
cuts is necessary in order to asses the final accuracy with which will be measured. 

A different proposal has been to simply discard all events with five or more jets. 
It is clear that this can only be a partial solution at best. It excludes hard, non- 
collinear gluon emission but it does not suppress soft and/or collinear radiation whose 
effects need in any case to be studied. Furthermore such a procedure would introduce 
a dependence of the mass measurement on the invariant mass cut which separates 
four-jet events from the events with a larger jet multiplicity. 

QCD radiation can modify the characteristics of the final state like the distributions 
of jet momenta, jet-jet invariant masses and opening angles from which the value of the 
W mass will be extracted. These semi-inclusive quantities, on the other hand, are of 
calorimetric type, that is infrared and collinear safe. The measured invariant masses, as 
an example, do not change if a massless parton of momentum p splits into two partons 
of momentum Xp and (1 — X)p or if an additional soft gluon is radiated by any of the 
hard partons. As a consequence they can be reliably studied in perturbative QCD. It 
is well known that differential distributions, particularly when the the phase space for 
gluon emission is restricted by experimental cuts, can be more sensitive to higher order 
corrections than fully inclusive quantities like total cross-sections in which virtual and 
real contributions cancel to a large degree. It is therefore necessary to include higher 
order QCD effects into the predictions for WW production and decay in a way which 
allows to impose realistic cuts on the structure of the observed events. In this paper 
we present the calculation of the complete 0{af) corrections to the three diagrams 
which describe the basic process e~^e~ — >• WW — > qiqiq2<i2- We will discuss only fully 
hadronic decays but the extension to semileptonic decays is trivial. 

Recently it has been pointed out §|, 05 fhat cross-talk between the decay of 
the two W^'s might take place also at energy scales much smaller than those typical of 
jets through color reconnection phenomena and Bose-Einstein correlations. This could 
weaken the notion of two separate decays and cast doubts on our ability to reconstruct 
the masses of the original sources from the decay products, producing potentially large 
uncertainties. Our understanding of these issues is, at present, extremely poor, and 
the predictions are strongly model-dependent. In this letter we will leave aside these 
non-perturbative issues. For a discussion of some aspect of colour reconnection from a 



perturbative point of view see & 10 



Calculation 

In n dimensions, = 4 — 2e, using dimensional regularization, the tree-level amplitude 
squared for the decay W ^ qq is: 

\Mo{qq)\^ =AKs{l-e) (3) 

with K = Nce^Cg and, labelling the momenta with the particle names, s = {q + qf. 
The corresponding amplitude squared for the decay W qqg is: 

|Mte)r = 8/r(l-e){(^-H£-Hi + £5 + il)-.(^H + il + 2)| (4) 

[\SiS2 Si S2 Si S2J \Si S2 J) 

with K' = KCpgslJ^'^^ s = {q + q + gf, Si = {q + gf and S2 = {q + gf. 



The phase space can be expressed as follows: 

dLips{Q q + q + g) = dLips{Q —^q + q)x 

s^'-^ {Air 



167r2 rd - 



dSidS2{{s - Si - S2)sSiS2} " (5) 



after integrating over the azimuthal angle in the center of mass of q and g with q on 
the 2;-axis. 

Integrating over Si and S2 in the region S where min(si, S2) < A one obtains: 

\M^{qq)f = fj^2 rlr, e) / £ ^"^^'^ {{s - s, - 82)33,32}-' \Miqqg)f (6) 

- mqq)f C.i-2^ {hi-^^'v^^n (^) (2, - 1,^ - ^) (7) 

+l + lv-2U2{v)-lu2{l)] 

= \Mo{qq) I' Cp^^i I ^ + A _ ^ - ^ In + I - ^ Li2(l) } + 0(r^ In r/) (8) 

where rj = A/3 and ^ = (^^^ e-e{7-in(47r))_ Li2(x) is the standard dilogarithm with 
Li2(l) = vr^/G. A separates the region of soft and collinear emission from the region 
where gluon radiation is considered hard. The dependence of the result on rj will be 
discussed later in detail. 

The one-loop contribution to the W decay is 

Mviqq) = Mo{qq)Cp£-^^ " I " ^ + ^ ^''^^^ " + D } ' 

Combining the virtual contribution (|^) with the integral of the real-emission cross sec- 
tion over the soft and collinear region the following simple expression is obtained 
for the W decay at C(a/): 

\M{W^qq)f = |Mo(gg)|V2Re(M*(gg)Mv^(gg)) + |M2(gg)|'+|ik%(gg)|2(10) 
= \Mo{qq)f{l + F{r])) + \M^{qq)f (11) 
= \Mo{qq)f{l + F'{7^)) + \M^{qq)\' + 0{r]lnr]) (12) 



where 



CpOis i 2 , ^ f V \ fn 1 2 3 



F{v) = ^ |- In^ + In ^-^j (2^ - -v' - -) (13) 



i + ^^_2Li2(r/) + Li2(l)}, 



Cf^s f ,_2 _ 3 , 1 



= -^[-ln'r^--\n7i-- + U2il)y (14) 

In the previous formulae \M^{qq)f is defined in analogy to (^ as the integral of 
the real-emission amplitude squared \M{qqg)f over the hard-gluon region E where 



min(si, S2) > A. The expression for the sum of the virtual and soft-colhnear contribu- 
tions has been derived many times before in various contests (see for example |]12, 



and we are in agreement with previous results, possibly after a simple coupling redefi- 
nition. 

Exploiting the universality of soft and collinear divergencies these results can be 
directly carried over to e"^e~ — > WW — > gigig2^2- At 0{af) there is no interference 
between the decays of the two W's and the amplitude for e~^e' WW Q'i^iQ'2^2fi' 
splits into two orthogonal terms Mi and M2 describing the emission from the qiqi 
pair and the ^2^2 pair respectively. The quantity s in ([10| ) in the two terms has to be 
identified with the invariant mass sig of the qiqig system and the invariant mass S2g of 
the ^2^25' system respectively. Choosing rji = Ai/sig = rj2 = /S.2/s2g = t] we can write: 

\M{WW^Aq)\'^ = \Mof + 2Re{M*Mv) + \Mj^f + \M^\^ (15) 
= \Mof{l + 2F'{r])) + \M^f + 0{7]\nr]) (16) 

In (|T5|) Mq is the tree-level amplitude and My the one-loop C(a/) amplitude for 
e~^e~ WW Aq. {M-^f and are obtained integrating the tree-level amplitude 

squared describing e~^e~ —>■ WW — ^ qiqiq2<i2g over the gluon variables as in (||). |Ms|^ 
is equal to the sum of the integral of |Mi|^ over the region Si where min((gi -l-^f)^, (gi -|- 
gY) < Ai and of the integral of |M2|^ over the region E2 where min((g2+5')^, (92+5^)^) < 
A2. IM^I is equal to the integrals of the square of Mi and M2 over the complementary 
regions. These latter integrals and the tree-level matrix element need to be calculated 
only in four dimensions. 

The real emission matrix element squared \M[X + g)\^ factorizes into the product 
of the corresponding emissionless expression |M(X)|^ times a universal function which 
describe the radiation of a gluon only in the soft and collinear limit. Therefore in the 
expression for the sum of the virtual and soft-collinear corrections there are unavoidable 
inaccuracies of order 77, related to the approximation of neglecting the gluon energy or 
emission angle inside |M(X)|^, and the use of the approximate expression (0) produces 
no loss in precision. In addition it must be remembered that we are effectively lumping 
the soft-collinear region in the phase space point corresponding to the absence of gluon 
radiation. Therefore, in order to obtain an accurate determination of the total cross 
section and in order to avoid disrupting the shape of the distributions of the different 
observable quantities a small value of 77 must be chosen. In all our figures we have used 
?7 = 1. X 10~^ and we have checked that our result depend extremely weakly on the 
value of 7] in the range 5. x 10^^ < < 1- x 10^^. 

Taking advantage of the fact that close to the appropriate soft-collinear regions, 
which contribute most to the result, the singular part of the integrands behaves like 
\M{qqg)\^ (|) we can use y\^i = ln(gj -|- g)"^ and y2,i = ln(gj + g)'^ as integration vari- 
ables. In this way it becomes possible to evaluate the real emission cross section using 
a reasonable amount of CPU time, producing theoretical predictions of an accuracy 
adequate to the expected precision of the measurement of the W mass. 

In the case of real emission, in order to extract two candidate masses from the event. 



The formalism is discussed in detail in ref. Il] 



the two partons with smallest 



2E,E,(l-cos%) 
^ s ^ ^ 

are merged, summing the corresponding four-momenta. This corresponds, at the 
present order in perturbation theory, to forcing all events to four clusters before at- 
tempting to reconstruct the masses of the two W bosons. In this paper we have re- 
stricted our attention to the simple JADE scheme and it remains to be studied whether 
a different reconstruction scheme [ll5| , p!6|] might prove more appropriate. Finally, the 



reconstructed masses are determined, out of the three possible combinations, choosing 
the two two-jet pairs whose masses M^i and MFt2 minimize 

Am = \Mr, - Mw\ + \Mr2 - Mw\ . (18) 

It is obvious that the precise form of the metric which defines the distance in jet -space 
is to some extent arbitrary. For instance one could use 

A^, = {Mm - Mwf + (M^2 - Mwf ■ (19) 

We have checked that at ^/s = 160 GeV and at ^/s = 175 GeV the cross sections 
obtained with ([T8|) and those obtained with ([19D differ by less than 1%. 

The expressions presented in this section can be applied without modification to all 
four-quark neutral current processes e~^e~ ZZ, Z'y,'j'j qiqiq2<l2- 

In order to expose the effect of 0{as) corrections we have made a number of sim- 
plifying assumptions. (1) Initial-state-radiation (ISR) effects have been neglected. It 
would however be straightforward to include them using standard techniques |T^. (2) 



Coulomb corrections have not been included although, since they can be expressed as 
a multiplicative factor times the lowest order (LO) matrix element squared, their in- 
clusion would be rather simple |]18[. (3) Quark masses have been neglected since the 



contribution of 6-quarks is severely suppressed by the smallness of the Vbc element of 
the CKM matrix and the c-quark mass is so small compared with Mw- (4) We have 
not taken into account non-resonant electroweak contributions and QCD four-jet back- 
grounds |19|, ^. The former has been shown to be well below the per-cent levelQ 
(with the possible exception of final states including electrons or electron-neutrinos). 
The latter, which is dominated by qqgg production, is strongly suppressed when two 
jet pairs with masses close to M]y are required. In ref. the QCD background, with 
Vcut = 0.01 and requiring two jet-jet masses within 10 GeV of M^y, has been estimated 
at about 0.1 pb for y/s = 161 GeV. Since the QCD background scales only as this 
estimate can be considered valid, as a first approximation, in the full range of operation 
of Lep II. However this procedure for reducing the four-jet background has far reaching 
consequences on the WW signal as will be discussed at length in the following. (5) 
Non-perturbative colour-reconnection and Bose-Einstein effects have not been consid- 
ered. (6) We have ignored all experimental uncertainties in the reconstruction of jet 
momenta and directions which will have to be studied with a full detector simulation. 



■^However they are intimately related to the gauge invariance of the theory. For a detailed discussion 
of this issue see ^ p2[ . 



All tree-level matrix elements have been computed, in the unitary gauge, using the 
formalism presented in ref. pSf with the help of a set of routines, called PHACT ^4 



which generate the building blocks of the helicity amplitudes semi-automatically. The 
matrix element for e~^e~ —>■ WW — > qiqiq2(j2g has been first computed in ref. [^. Our 
result has been checked against the corresponding amphtude generated by Madgraph 



25|. 

The numerical values of the Standard Model parameters used in the calculation are 
given in Table II. 



Results 

Our results are presented in fig.l, fig.2. Table III and Table IV. The total hadronic 
cross sections at tree-level and at next-to-leading-order (NLO) in QCD are given 
in the first column of Table III for ^/s = 160 GeV and 175 GeV. These results are 
sensitive to events with jet pairs of arbitrarily large and small invariant masses within 
the kinematical limits. In practice events with jet-jet masses too far from the W mass 
would not be accepted as WW events. Cuts of this sort have been suggested in order 
to eliminate most of the four-jet QCD background. We therefore require that each of 
the reconstructed masses, selected according to (p!8D, satisfies 

\MRi-Mw\<S, 2 = 1,2. (20) 

In the second and third column of Table III we show the cross section obtained with 
(5 = 30 GeV and 6 = 10 GeV respectively. In parentheses we give the cross sections 
obtained when all events which are identified as containing five jets, with ycut = 1. x 10^^ 
in the JADE scheme, are discarded. The ratio of the corresponding five-jet cross section 
to the total hadronic cross section 

_ a{WW^5j) 
^ a{WW had) 

is about 22% at y/s = 175 GeV with no mass cut. 

Let us mention that in the so called ADLO-TH set of cuts, the semi-realistic set 
which the Working Groups on LEP II Physics have agreed on in order to provide a 
common ground for the comparison between theoretical calculations and the simulations 
performed by the experiments, two jets are considered as distinguishable if the invariant 
mass of the pair is larger than 5 GeV. This corresponds to ycut ~ 1. x 10^'^. However 
at this value of ycut the five-jet cross section is larger than the total hadronic cross 
section and, as a consequence, the cross section for events with at most four jets is 
negative. This can be readily seen from eq. (0). If we neglect for simplicity the 
interplay between the two decays, we have that ycut = 1. x 10~^, since ^/s ~ 2M\y, 
corresponds to r/ ~ 4. x 10"'^ and one finds that, with = .117, F'(0.004) < —1. 

In the narrow width, or stable W^s, approximation the QCD correction to the total 
cross section is simply twice the correction to the W decay, namely the Born cross 
section gets multiplied by (1 + 2as/7r). In principle, for off-mass-shell W^s, one could 
expect additional terms of order {asTw)/{T!'Mw) ~ 1. x 10~^. In order to study effects 
of this order of magnitude, all matrix elements have been integrated using VEGAS with 



(21) 



a relative precision not worse than 1. x 10~^. We have found that the corrections to the 
total hadronic cross section for unstable W^s do not differ from those for stable W's by 
more than one part in a thousand. Therefore the narrow width approximation appears 
to be perfectly adequate for the QCD correction to the total cross section 0. This result 
however relies on the delicate cancellation between real and virtual contributions, and 
is strongly modified by any constraint which decreases the phase space available for real 
emission. This is evident from the second and third column of Table III. At ^/s = 175 
GeV , with 6 = 10 GeV, NLO QCD corrections decrease the corresponding tree-level 
cross section by about 33%, which is to be contrasted with the increase of about 7% they 
produce in the total hadronic cross section in column one. Even with the much milder 
cut 6 = 30 GeV the LO result is decreased by about 11%. At = 160 GeV, about 
200 MeV above the nominal threshold with the value of Mw adopted in this paper, 
the NLO result with S = 10 (30) GeV is 45% (15%) smaller than the LO prediction 
without mass cuts. The corresponding cross section is decreased by approximately 0.8 
(0.3) pb. It should be remembered that a change in the total cross section of 0.2 pb 
corresponds to a shift of about 100 MeV in the measured W mass ||^. Therefore QCD 
corrections, in contradiction to the naive expectations based on the behaviour of the 
total cross section, are potentially as large as the corrections due to ISR and much 
larger than those produced by the Coulomb interaction, even though they only affect 
about half of the total e^e~ W^W~ cross section. Obviously, any modification 
of the LO cross section which can be reliably computed can be incorporated into the 
theoretical predictions which are fitted to the data and do not degrade the expected 
accuracy of the measurement of Myy. They however directly affect the statistical error 
and the uncertainty in the theoretical analysis has to be included in the evaluation of 
the systematic error. 

The behaviour of the LO and NLO hadronic cross section as a function of the center- 
of-mass energy can be found in fig.l. As in Table III we present the total hadronic cross 
sections at tree-level (dashed curve) and at NLO (continuous curve). The long-dash- 
short-dash curve and the dash-dot curve show the NLO cross sections obtained with 
6 = 30 GeV and 6 = 10 GeV respectively. The dotted line gives the tree-level cross 
section with 6 = 10 GeV. The analogous cross section with 6 = 30 GeV, as expected 
from the fact that in this case 6 ^ Tyy, differs very little from the total cross section 
and we omit it. 

The effect of the mass cut (^0|) reaches a maximum at about ^/s = 155 GeV and 
decreases with increasing center-of-mass energy. We notice that at ^/s = 150 GeV, 
with 6 = 10 GeV, the NLO result becomes larger than the LO one. This is due to 
the fact that typically, below threshold, one of the W^s is on mass-shell while the 
other is highly virtual. Therefore at tree-level most of the events fail the cut. In 
this case the rearrangement of invariant masses introduced by QCD radiation blurs the 
separation between on-mass-shell and off-mass-shell W^s and reduces the effectiveness 
of the mass cut. From fig.l we read that the NLO order result with 5 = 30 GeV is 
always considerably smaller the tree-level prediction, confirming that even moderate 
restrictions on the structure of the final state can produce sizable effects. Choosing 
S = 10 GeV decreases the available event rate by a factor close to 50% and might 

^As a consequence, the QCD corrected cross section for W^W^ — > qqiiy, if no restriction is imposed 
on the hadronic part of the final state, is cro(l + as/n) where ctq is the Born result. 



result in a substantial increase in the statistical error of the VT-mass measurement. 
Just above threshold the QCD corrections depend only weakly on the center-of-mass 
energy. At = 160 GeV they amount to 45% (15%) with 6 = 10 (30) GeV while at 
^/s = 165 they become 38% (14%). Therefore we do not expect the determination of 
the optimal energy for the threshold measurement of Mw of ref. to be modified by 
QCD effects. 

In fig. 2 we present the distribution of the average of the two reconstructed W masses 

M=^{Mni + Mn2). (22) 

which has been extensively used as an estimator of M\y. At tree-level the difference 
between M and the average of the masses of the two virtual W^s, namely the masses 
which one would reconstruct if the quarks could always be correctly paired, is extremely 
small. Clearly one could consider more sophisticated approaches, for instance using the 
mass distribution in the {Mjn, M^2) plane but this analysis is beyond the scope of this 
letter. The two upper curves in fig. 2a are obtained when all events are retained, even 
when one or both the candidate masses are far from Mw- The results for 6 = 30 GeV 
and 6 = 10 GeV are shown in fig. 2b and fig. 2c respectively. The dashed lines are the 
LO result while the NLO predictions are given by the continuous lines. The NLO cross 
section for events with at most four jets, at i/cut = 1. x 10~^ is shown by the dotted lines. 
When no cut on |M^j — M]v\ is imposed, at NLO the peak cross section is strongly 
reduced with respect to the LO distribution and a large tail at low masses is generated. 
As S decreases the tail gradually disappears while the number of events contained in 
the peak at about the W mass decreases only slightly. It is precisely the distortion 
of the average-mass tree-level spectrum that explains the dramatic effect of the mass 
cut ( PPD at NLO. We see that the low-mass tail persists also in the lower multiplicity 
sample and therefore it is not predominantly populated by five-jet events. 

In an attempt to quantify the modifications of the average mass distribution which 
are induced by QCD corrections we present in Table IV the average deviation, (M — 

Mw) and the root-mean-square (RMS) deviation, {(^M — Mw^ of the mass dis- 
tribution with respect to Mw in GeV. In each column the left (right) hand side value 
refers to the LO (NLO). In the first row all events are accepted while in the second and 
third row the two jet-jet masses satisfy jM^j — Mw\ < S,i = 1,2 with 5 = 30 GeV and 
S = 10 GeV respectively. The results in the last three rows are obtained discarding all 
events which are recognized as five jets with i/cut = 1. x 10~^ in the JADE scheme. 

Table IV shows that already at tree-level the average reconstructed mass M is 
smaller than Mw- The difference AM = Mw — M depends on the mass cut. AM can 
be larger than 300 MeV and decreases with 6 down to about 40 MeV at 5 = 10 GeV. 
This behaviour can be observed directly in fig. 2. The asymmetry which is evident in 
fig. 2a decreases as the mass cut is tightened. When no mass cut is imposed the mass 
distribution is so deeply altered that it is almost meaningless to compare the tree-level 
two lowest order moments with those at NLO. This is still partially true for 6 = 30 
GeV. The low mass tail is non-negligible and considerably increase the mass shift from 
.274 MeV to 1.70 GeV and the RMS width of the distribution, from 2.39 GeV to 5.11 
GeV. When the mass cut is 5 = 10 GeV AM increases from 41 MeV to 125 MeV while 
the RMS deviation increase from about 1.7 GeV to slightly more than 2 GeV. 



We believe that the impact of the results presented in this paper on the precision 
measurement of M\v^ both at threshold and through direct reconstruction, should be 
carefully assessed taking into account the other known sources of large corrections, 
namely ISR and Coulomb interaction, and the experimental uncertainties. In our 
opinion an improvement of our understanding of the four-jet QCD background and 
a refinement of the strategies to reduce it are required . As for many other effects, if 
the theoretical analysis can be made sufficiently precise and, at the same time, accurate 
data are available, through a careful extrapolation of LEP 1 results, or can be obtained, 
possibly through a dedicated run below threshold, QCD backgrounds could become part 
of the predictions without compromising the accuracy with which the 1^-mass will be 
measured. Further work will be necessary to determine which combination of cuts will 
result in the highest accuracy on M^, balancing the smaller numer of events with a 
decreased level of background. 

Conclusions 

We have presented the calculation of the complete 0{aj) corrections to e+e" WW 
?i^i?2^2 in a form which allows to impose realistic cuts on the structure of the observed 
events. It has been shown that QCD radiation modifies in an important way the 
distribution of jet-jet invariant masses and therefore that it strongly influences the 
measurement of the W mass from direct reconstruction of the decay products. For 
total hadronic cross-sections virtual and real contributions cancel to a large degree if 
all events are retained. If however only events with two jet-jet masses close to Mw are 
accepted, higher order QCD corrections can substantially decrease the cross section. 

A number of issues have not been considered in this paper and will have to be 
addressed in order to obtain a more complete theoretical understanding of the pro- 
duction of WW pairs at Lep II and to determine the best strategies for the precision 
measurement of the W mass. (1) The precise choice of the distribution or combination 
of variables from which to extract Mw- (2) The influence of the jet reconstruction al- 
gorithm including hadronization. (3) Initial state radiation. (4) Coulomb corrections. 
(5) Non resonant background effects. (6) Colour reconnection effects. As already men- 
tioned, initial state radiation and Coulomb corrections are relatively straightforward 
and an improved analysis which include both effects is well under way. 
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Figure Captions 



Fig.l Hadronic cross section for e^e~ — > l^"*"!^" as a function of ^/s. The dashed hne 
and the dotted hne are tree-level predictions. For the latter we have required 
\MRi — Mw\ < 10 GeV, i = 1, 2. The continuous line is the total hadronic NLO 
result. The long-dash-short-dash curve and the dash-dot curve show the NLO 
cross sections obtained with 5 = 30 GeV and 5 = 10 GeV respectively. ISR and 
Coulomb effects are not included. The parameter values are given in Table II. 

Fig.2 Distribution of the average of the two reconstructed mass M = ^ {Mm + Mfa)- 
The dashed lines are the tree-level predictions. The continuous lines and the 
dotted lines include NLO QCD corrections. For the latter five-jet events, defined 
using Ucut = 1. X 10~^ in the JADE scheme, have been excluded. In fig.2a no 
cut on iM/jj — Mw\ has been imposed. In fig. 2b and fig. 2c we have required 
|Mffi - Mw\ < 10 GeV i = 1, 2 with 5 = 30 GeV and 5 = 10 GeV respectively. 
ISR and Goulomb effects are not included. The parameter values are given in 
Table II. 



Table Captions 

Table I Branching ratios of the WW pairs at tree-level. The decay channels involving 
r's are particularly challenging and are usually analized separately. 

Table II Parameters used in the numerical part of the paper. 

Table III Total hadronic cross sections. For each energy the upper line refers to the 
LO and the lower one to the results at NLO in QCD. In the first column all 
events are accepted while in the second and third the two jet-jet masses satisfy 
|Mffi - Mw\ < = 1,2 with 5 = 30 GcV and 5 = 10 GeV respectively. The 
numbers in parentheses are obtained discarding all events which are recognized 
as five jets with ycut = 1. x 10~^ in the JADE scheme. 

Table IV Average deviation and RMS deviation from My/ in GcV. In each column 
the left hand side value refers to the LO and the right hand side one to the results 
at NLO in QCD. In the first row all events arc accepted while in the second and 
third row the two jet-jet masses satisfy \MRi — Mw\ < 5,i = 1,2 with 5 = 30 
GeV and 5 = 10 GeV respectively. The results in the last three rows are obtained 
discarding all events which are recognized as five jets with ycut = 1. x 10~^ in the 
JADE scheme. 
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Table I 



parameter 


value 


Mz 


91.173 GeV 


Mw = Mz cos 6*^^ 


79.905 GeV 


Tz 


2.49 GeV 




2.08 GeV 




128. 


sin^ 


0.2319 




.117 




.38937966 10^ p^GeV^ 



Table II 







5 = 30 GeV 


5 = 10 GeV 


= 175 GeV 
{Vcut = 10-') 


6.749 


6.725 


6.393 


7.251 (5.638) 


6.169 (4.718) 


4.277 (3.383) 


= 160 GeV 


1.784 


1.757 


1.513 


1.916 


1.513 


0.980 



Table III 



x/i = 175 GeV 


(M - Mw) 


((m-m^)V 




LO NLO 


LO NLO 


no cut 


-0.340 -4.82 


2.645 10.65 


5 = 30 GeV 


-0.274 -1.70 


2.390 5.11 


5^10 GeV 


-0.041 -0.125 


1.679 2.02 


ycut = 0.01 


-0.340 -3.88 


2.645 9.88 


S = 30 GeV, y^t = 0.01 


-0.274 -1.09 


2.390 4.15 


5 = 10 GeV, Vent = 0.01 


-0.041 -0.092 


1.679 1.78 



Table IV 
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